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Abstract 

Wc consider the Kuramoto model of globally coupled phase oscillators with time-delayed inter- 
actions, that is subject to the Ornstein-Uhlenbeck (Gaussian) colored or the non-Gaussian colored 
noise. We investigate numerically the interplay between the influences of the finite correlation time 
of noise r and the time delay on the onset of the synchronization process. Both cases for identical 
and nonidentical oscillators had been considered. Among the obtained results for identical oscilla- 
tors is a large increase of the synchronization threshold as a function of time delay for the colored 
non-Gaussian noise compared to the case of the colored Gaussian noise at low noise correlation time 
r. However, the difference reduces remarkably for large noise correlation times. For the case of 
nonidentical oscillators, the incoherent state may become unstable around the maximum value of 
the threshold (as a function of time delay) even at lower coupling strength values in the presence of 
colored noise as compared to the noiseless case. We had studied the dependence of the critical value 
of the coupling strength (the threshold of synchronization) on given parameters of the stochastic 
Kuramoto model in great details and presented results for possible cases of colored Gaussian and 
non-Gaussian noises. 

PACS numbers: 05.40.-a, 05.45.Xt, 89.75.-k, 02.30.Ks 



I. INTRODUCTION 



Since the pioneering works on coupled phase oscillators by Winfree jlj] and Kuramoto 
syuchrouization ■„ „o,J,„ear systems has bee,. syste,„at,ca..y studied and attracted much 

|. The systems where 

the phenomena of synchronization 



attention. Recent reviews on the developments can be found in [3|-_^ 

lave been observed include biological clocks chemical 
7H9| , coupled random frequency oscillators |lO[ , cardiores- 



oscillators [2| , coupled map 



attices 



piratory coupled system etc. 

The model first introduced by Kuramoto is one of the basic models that describes the 
synchronization process when initially independent oscillators begin to move coherently. It 
was thoroughly studied and successfully applied in several systems which were modeled by an 
ensemble of coupled phase oscillators |6|. Another model of interest is the Kuramoto model 



with a tim^delay |lM5|. The mode, shows a number of interesting phenomena ineh.dirrg. 

e.g., the effect of bistability as discovered in [1^ where the Kuramoto model with time-delayed 
interactions was also considered to be subject to a white noise. Note that both time-delayed 
interactions and noise play very important role in nature (see, e.g., fl6|). 

In this paper we consider the Kuramoto model of globally coupled (all-to-all) phase oscil- 
lators with a time delay and subject to the Ornstein-Uhlenbeck (OU) colored noise [iTi that 
is a Gaussian process with a finite correlation time, and a non-Gaussian colored noise 18N21|. 
The focus in our study is on the interplay between the influence of noise and the time delay 
on the synchronization process. Our previous study on the stochastic Kuramoto model 21 | 
showed that the influence of the OU noise qualitatively differs from the case of white noise as 
the former allows for the full synchronization despite the fact that the system is subject to a 
noise, see e.g. Fig. 3 in [21]. 

Generally speaking, the time delay introduces de-phasing among the oscillators. This de- 
phasing is to interfere with the intrinsic correlations caused by the finite correlation time 
of the noise. As a result, we have found that the de-phasing plays important role in the 
dynamics of the system when the intrinsic correlations are small. We also investigate the 
effect of the time delay on the onset of synchronization for different noise strengths for both 
OU and non-Gaussian noises. But before doing that we show that the effect of time delay for 
the noiseless Kuramoto model is qualitatively equivalent to the effect of frequency fluctuations 
of the phase oscillators. Meanwhile we demonstrate that for the stochastic Kuramoto model 
the critical coupling grows nonlinearly with the increase of the time delay that is due to the 
additional de-phasing caused by the time delay. We also investigate how the transition from 



the Gaussian to a non-Gaussian noise changes the dynamical properties of the system. 



II. STOCHASTIC KURAMOTO MODEL WITH TIME-DELAYED INTERAC- 
TIONS 

Let us consider the stochastic Kuramoto model with time-delayed interactions. This model 
describes coupled phase oscillators with dynamics governed by the following equations 

^ = + ^ E - ^'^) - ^^(^)] + ^^('^)' (1) 

where 9i and Ui are, respectively, the phase and the frequency of the i-th oscillators {i = 
1, . . . , A^), e is the coupling constant, and is the time delay. The independent noise processes 
rjiit) are governed by 

dr]i{t) ^ d 

—1— = — + — im- (2) 

dt T drji T 

The potential function is 

(^^) = V{^T) ^""^^ + " ^^''^'/^^ 

with a = t/D. ^i(t) is the Gaussian white noise process defined via 

mMt')) = 25,j5it-t') 

and (^i(t)) = 0. D and r measure the intensity and the correlation time of the noise process. 
Shortly we will discuss more about them. In the noise term in Eq. (1) it is apparent that in 
the present problem we assume the homogenous diffusion of phase oscillators. 

The form of the noise rji allows us to control the deviation from the Gaussian behavior by 
changing a single parameter p. For p = 1, Eq.(l2]) becomes 



(3, 



which is a well-known time evolution equation for the OU noise process 17|] for which the 
auto correlation function is given by 

< mitMt') >= - exp(-^^) . (4) 
r T 

Thus D and r are the noise strength and correlation time of the OU noise. The factor 
[1 + a{p — l)rif/2] in the first term of Eq.(l2]) leads to producing colored non Gaussian noise 



whose effective noise strengtli and correlation time are different from D and r. However, to 
make the present paper self-consistent we would like to mention here salient features of the 



non- Gaussian noise only ISj. 



The stationary probability distribution of the noise process is given by 

1 



p-1 

(5) 



where Zip is the normalization constant which equals to 

.1 — ^ ^ 



-1 



ri(i/(p-i)-i/2) 

(6) 



a{p-l) ri(l/(p-l)) 

with Fi being the Gamma function. This distribution can be normalized only for p < 3. Since 
P{rii) is an even function of r^j, the first moment, (t^j), is always 0, and the second moment 
(variance) is given by 

^^-^ = 7(5^ ' 

which is finite only for p < 5/3. Furthermore, for p < 1, the distribution has a cut-off and it 
is only defined for \'qi\ < r]ic = \J~^^^- 

However, it is difficult to determine the auto correlation function (ACF) of the non Gaus- 
sian noise r/, exactly. To have an idea about this we have calculated it numerically and 
presented the result in Fig. 1. The two time correlation function for the non Gaussian noise 
(solid curve) is fitted well by bi-exponentially decaying function < 'r]i{t)r]i(0) >= ae~*/'^^ + 6e^* 
having two correlation times 37 and 1, respectively, for r = 1.0, where a and b are constants 
related to the noise intensity. In the same figure we have plotted numerically calculated ACF 
for Gaussian colored noise. It exactly mimics the function given in Eq.(jl]). Figure 1 clearly 
shows that the effective noise strength and correlation time for non Gaussian noise are greater 
compared to those of colored Gaussian one. The correlation time of non-Gaussian noise Tp at 
';he stationary regime of the process rii{t) diverges near p = 5/3 and it can be approximated 



18| over the range l<p<5/3as 



rp~2r/(5-3p) . (8) 

This equation is qualitatively consistent with Fig.l. However, for this approximate corre- 
lation time, Eq.(I7]) becomes 



Equations (I7}|9]) imply that when p — )■ 1, we recover the hmit that T]i is a Gaussian colored 
noise since at this limit Eqs.f l7f9|) correspond to the variance of Colored Gaussian noise as 
given in Eq.(j4]) and Tp in Eq.(8) becomes r which is the correlation time of the Gaussian noise. 
Another check in this context can be obtain in the following way. In this limit the term in 
the square bracket of Eq. can be written as 

l + a{p- l)r,f/2 = exp(a(p - l)rfj2) . (10) 

Then Eq.Q becomes 

p(^) = i_exp(-«r/2/2) , (11) 

with Zi = A/vr7a, which is a Gaussian distribution function. However, Eq.([7]) shows that 
for a given external noise strength D and the noise correlation time r, the variance of the 
non-Gaussian noise is higher than that of the Gaussian noise for p > 1, i.e. {rjf^) > {rjf^). 
Similarly, Eq.(l8]) implies that Tp > t for p > 1. These are consistent with the message of 
Fig. 1. Before leaving this issue we like to emphasize that in the present study we have 
considered continuous distribution of the non-Gaussian noise which is more relevant to the 
natural systems rather than two-state or discrete distributions as mostly used in the literature 



22|. 



to study the noise driven dynamical systems 
The quantity of interest in the present study is 

1 ^ 

i=l 

which is the order parameter that measures the extent of synchronization in the system of 
phase oscillators. Its absolute value F determines the degree of synchronization. It can 
be seen that in case of all the oscillators having the same phase the quantity equals to one 
(F = 1) that corresponds to the full synchronization. The degree of synchronization is equal 
to zero (F = 0) when all the oscillators are independent and have different phases. B defines 
an average phase of the oscillators. For the initial distribution of frequencies we choose the 
Lorentzian distribution 



III. RESULTS AND DISCUSSION 



We have computed the order parameter T as well as the critical coupling strength numer- 
ically since it is very difficult to study the problem analytically due to the colored noise and 
nonlinearity in terms of rji in its time evolution equation. Using Heun's method (a stochastic 
version of the Euler method which reduces to the second order Runge-Kutta method in the 
absence of noise) we have solved the Eqs. ([1]) and ([2]) simultaneously [23|. Based on the above 
method, we have studied time evolution of = 5000 coupled phase oscillators. In Ref. |21| . 
we have showed that such a value is large enough so that the calculated quantities can 
represent those for systems in the thermodynamic limit. 

For a given frequency distribution of the phase oscillators and environment of specific 
characteristic there exists a threshold value of the coupling strength (ec) above which the 
coherent (synchronized) state is the stable state. If the e is smaller than ec, then the incoherent 
state is the stable one. ec is called critical coupling strength. To calculate the value of ec, we 
ffist obtain the dependence of ensemble average of F (< F >) on e. From this we determine 
numerically the derivative of < F > with respect to e at stationary state. The e corresponding 
to the maximum value of the derivative is identified as the critical coupling strength. Our 
numerical study satisfies all the known limiting results, i.e., (i) ec = 2 A in the absence of noise 
{D = 0) and = (ii) ec = 2(A + D) for white noise and = (results are not shown here). 
We would also like to mention here that we chose stationary time to be a linear function of 
time delay (r^) with the proportionality constant 100 and integration step length h = 0.01. 

We start our numerical study to determine how stability of both coherent and incoherent 
states depends on the time delay. Stability of these states of white-noise driven coupled 
oscillators in the Kuramoto model with a time delay has been studied by Young and Strogatz 
[l^ . They have considered two cases. In the ffist case, all the oscillators are identical. 
In the second case, they have chosen Lorentzian distribution, Eq.f llSp . of frequency of the 
coupled oscillators. In the present paper, we extend this study to colored-noise driven coupled 
oscillators. The noise may be OU Gaussian or non Gaussian in characteristics. However, the 
critical coupling strength as mentioned above indicates that up to that value of coupling 
strength, the incoherent state can survive at long time. This range corresponds the stability 
zone of the incoherent state of the coupled phase oscillators. Thus critical coupling strength 
is a measure to imply how stability of incoherent or coherent states of the coupled oscillators 
depends on time delay and other parameters of the system. To demonstrate the dependence 
of the critical coupling strength ec on the time delay for various noise properties, we have 



reproduced Fig.2 of Ref. IJ] for identical oscillators in the absence of noise (i.e. D = 0) 
and presented it by the solid curve in Fig.2. In the presence of noise the phase oscillators 
become non identical and the time delay is effective in the dynamics even at its low value 
compared to the case where noise is absent. Since the variance of the non-Gaussian noise 
is higher compared to the Gaussian one, the shift of the damped oscillating curve towards 
the larger critical coupling strength from the noiseless case, is very big for former than latter 
at low noise correlation time, r ^ 0.1. Their difference reduces remarkably at large noise 
correlation time, t ^ 2. The reason may be the following. Effective noise correlation time 
of the colored non-Gaussian noise is much higher than the Gaussian one at large r. Increase 
of noise correlation time leads to develop better phase relationship among the oscillators 
reducing the phase diffusion. Thus at large r the difference in critical coupling strength is 
small for a given time delay. It is apparent in Fig.2 that the interplay of noise correlation 
time and the time delay plays some constructive role to have a synchronized state particularly 
at large r and around the first maximum of the damped oscillation curve. It also shifts a 
little bit the position of the maximum towards the left and reduces the oscillation amplitude 
at large time delay values. Thus colored noise plays a role beyond the induction of diffusion 
behavior. 

Before going to the next case we would like to demonstrate the variation of critical coupling 
strength as a function of frequency of the coupled identical oscillators in the presence of time 
delay. In Fig. 3 we have presented this. It exhibits that the variation is periodic as a result 
of interplay of the time delay and the frequency of the oscillator. However, again it shows 
how noise correlation can reduce the stability of incoherent state. Because of higher effective 
noise strength of the non-Gaussian noise, the critical coupling strength for the case p = 1.5 
and r = 0.5 is comparable to white Gaussian noise. 

For the second case of non-identical phase oscillators, we have done calculations similar to 
those for Fig. 2 for identical oscillators, and the results are presented in Fig. 4. Again, the 
solid curve in Fig.4 is a reproduction of Fig. 4 of Ref.jl^. All the features of Fig. 2 in the 
presence of colored noise has appeared in Fig. 4. Here one important point to be noted is that 
around the maximum, the incoherent state is unstable even at lower coupling strength in the 
presence of colored noise compared to the case without noise. Thus colored noise can effect 
the nonlinear coupling among the phase oscillators. It is quite similar to the modification of 
the dynamics by the colored noise in the presence of nonlinear potential 2^, l25 |. 

In the next step we demonstrate how the critical coupling strength varies as a function 
of noise parameters. First, we consider the dependence of ec on the noise strength (D). To 



identify the signature of the time delay in this context we have determined a ratio of critical 
coupling strengths edrd = 0)/^c{Tdm) for a given noise intensity, where Tdm is the location of 
the first maximum in the axis. For identical oscillators we have chosen Tdm = 2 from Fig. 
2 and for nonidentical oscillators we have considered — 1 from Fig. 4. It is not difficult 
to anticipate that this choice will incorporate the maximum effect of time delay. However, 
the ratio had been calculated for different values of the noise strength for the white noise 
and is presented in Fig. 5. At first it rapidly increases both for identical (dotted curve) and 
nonidentical (solid curve) oscillators then the growth rate slows down. The initial growth 
rate is higher for identical oscillators compared to nonidentical oscillators because at low 
noise strength edra — 0) is close to zero for the former and it grows at a faster rate. If the 
noise strength is appreciably large then the ratio is higher for identical oscillators than for 
the nonidentical case even if Tdm is greater for the former than the latter. Thus time delay is 
more effective in the dynamics for nonidentical oscillators compared to the case of identical 
oscillators. Similar features are also observed for the colored noise and are demonstrated in 
insets (a) and (b) for the Gaussian and non-Gaussian noises. At large noise intensities the 
ratio of critical coupling strengths is smaller for the colored Gaussian noise compared to other 
noises for the case of nonidentical oscillators. Thus on switching from the white Gaussian to 
colored Gaussian or non-Gaussian characteristics of noise the time delay begins to play an 
important role in the dynamics. 

The variation of the ratio of critical coupling strengths with the non-Gaussian parameter 
p is demonstrated in Fig.6. It shows that the ratio increases as a function of p for both cases 
of identical and nonidentical oscillators. The rate of growth is higher for the former than the 
latter case. The ratio is lower for the identical oscillators compared to the case of nonidentical 
oscillators at low p. However at large p the situation turns around, it becomes inverse. These 
aspects can be understood if one keeps in mind that the effective noise strength increases as 
p grows. Therefore Fig. 5 is essentially similar to Fig.6 and it can be explained in the same 
way as we have explained the results presented in Fig.6. 

Finally, in Fig. 7 we have presented the variation of the ratio of critical coupling strength 
as a function of noise correlation time. It shows that the ratio decreases with increase of 
T. Since phase diffusion reduces as the noise becomes more colored, the critical coupling 
strength decreases as r grows both in the presence and the absence of time delay. The ratio 
decreases because in the absence of time delay the colored noise is more effective in the 
dynamics. Similarly, the rate of decrease is a little bit higher for identical oscillators than for 
nonidentical ones. These kind of features are also observed for the colored non-Gaussian noise 



and are presented in the inset of the Fig. 7. In both cases the ratio is higher for the identical 
oscillators compared to the nonidentical oscillators for a given noise correlation time. Thus 
it again supports the statement that the time delay is more effective in the dynamics for the 
latter than for the former. 

IV. CONCLUSION 

We have considered the stochastic Kuramoto model including time-delayed interaction 
with a delay time subject to both OU Gaussian or non-Gaussian colored noise with a 
correlated time r. The main focus of our study was on the interplay between the effects 
of finite correlation time r and the time delay and their influence on the onset of the 
synchronization. Our results can be summarized as follows. 

(i) The shift of the damped oscillating curve (for the critical couphng strength as a function 
of time delay r^) towards the larger critical coupling strength from the noiseless case is very 
large for the non-Gaussian noise than for the Gaussian one at low noise correlation time, 
r ~ 0.1. Their difference reduces remarkably at large noise correlation times, r ~ 2 (Fig. 2). 

(iii) Due to the present intrinsic correlations, the colored noise plays a role beyond the 
induction of diffusive behavior. 

(iii) The incoherent (unsynchronized) state may be unstable around the maximum even at 
lower coupling strength values in the presence of colored noise compared to the noiseless case 
(Fig. 4). 

(iv) The ratio of critical coupling strengths ec(rrf = 0)/ec{Tdm) first rapidly increases as 
a function of noise intensity D for both the identical and nonidentical oscillators and then 
slows down. The initial growth rate is higher for identical oscillators compared to the case 
of nonidentical oscillators. If the noise strength is appreciably large then the ratio is higher 
for identical oscillators than for the nonidentical case even if is greater for the former 
than the latter case. At large values of noise intensity the ratio of critical coupling strength 
is smaller for the colored Gaussian noise than for the other noises for the case of nonidentical 
oscillators (Fig. 5). 

(vi) The ratio of critical coupling strengths increases as a function of non-Gaussian param- 
eter p for both cases of identical and nonidentical oscillators. The rate of growth is higher 
for the former case than the latter one. The ratio is lower for identical oscillators compared 
to nonidentical oscillators at low p > 1. However at large values of p it becomes inverse (Fig. 
6). 



(v) The above ratio decreases with increase of r. The rate of decrease is higher for the 
case of identical oscillators compared to the one for nonidentical oscillators (Fig. 7). 

We anticipate that investigation of influence of time-delayed interactions and (generally 
non-Gaussian) noise in complex systems would lead to important insights into their stochastic 
dynamics, e.g., for the process of intercellular synchronization in biology 26(| (see also systems 
biology models mentioned as possible applications in \^^). The time-delayed interaction also 



27| 



plays an important role in a molecular model of biological evolution 
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FIG. 1: Two-time correlation function vs. t of both Gaussian {p = 1) and non-Gaussian (p = 
colored noises for parameters r = 1.0 and D = 0.5. 
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FIG. 2: Plot of the critical value of coupling strength vs the time delay for the identical 
oscillators for = 7r/2. 




FIG. 3: Plot of the critical value of coupling strengths (ec) vs the frequency of identical oscillators 
ojQ for D = 1.0 and = 2.0. 




FIG. 4: Plot of the critical value of coupling strength Cc vs the time delay for the nonidentical 
oscillators for the parameter set loq = 3.0 and A = 0.1. 




FIG. 5: Plot of the ratio of the critical value of coupling strengths edrd = 0)/ec{Tdm) vs the noise 
intensity D for the parameter set p = 1.0 and r = 0.0. For identical oscillators ujq = Tr/2 and that of 
nonidentical oscillators, ojo = 3.0, A = 0.5. Solid and dotted curves correspond to nonidentical and 
identical oscillators, respectively. The insets contain results for different parameter sets: (a) p = 1.0 
and r = 0.5, (b) p = 1.5 and r = 0.5 
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FIG. 6: Plot of the ratio of the critical value of coupling strengths edra = 0)/ec(Tdm) vs the non- 
Gaussian parameter p for the parameter set r = 0.5 and D = 0.25. For identical oscillators cjq = 7r/2 
and that of nonidentical oscillators ojq = 3.0 and A = 0.5. Solid and dotted curves correspond to 
nonidentical and identical oscillators, respectively. 
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FIG. 7: Plot of the ratio of the eritical value of coupling strengths {ec{Td = ^) / ^ciTdm)) vs the noise 
correlation time r for the parameter set p = 1.0 and D = 1.0. For identical oscillators ojq = 7r/2 
and that of nonidentical oscillators ujq = 3.0 and A = 0.5. Solid and dotted curves correspond to 
nonidentical and identical oscillators, respectively. In the inset, p = 1.5. 



